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1.  INTRODUCTION 


The  chemical  simplicity  of  hydrocyanic  acid  (HCN)  and  its  multimers  makes  them  attractive 
candidates  for  experimental  and  theoretical  study.  Due  to  its  extreme  toxicity,  HCN  must  be  handled  with 
caution.  Additionally,  HCN  is  reported  to  be  unstable  upon  storage  with  the  possibility  of  explosion  due 
to  polymerization,  presumably  in  forming  1,3,5-triazine  (C3N3H3),  also  known  as  sym-triazine 
(Migrdichian  1947).  To  explore  the  possibility  of  high-energy  release  in  the  formation  of  sym-triazine 
from  HCN,  more  information  is  needed  on  the  reaction  mechanism  and  barrier  height  to  formation. 
Unfortunately,  very  little  is  known  about  the  potential  energy  surface  (PES)  for  the  polymerization  of  HCN 
to  jym-triazine.  There  have  been  experimental  studies  on  the  reverse  reaction,  the  decomposition  of 
.^m-triazine  to  form  HCN  (Ondrey  and  Bersohn  1984;  Coates,  Chu,  and  Flynn  1984). 

sym-triazine  3  HCN  AH  =  43.2  kcal/mol  (I) 

Photodissociation  experiments  at  248  and  193  nm  showed  that  sym-triazine  decomposes  in  a  concerted 
manner  only  and  forms  three  HCN  molecules  (Ondrey  and  Bersohn  1984).  These  measurements  provided 
an  upper  bound  of  the  barrier  to  reaction  (I)  (1  IS  kcal/mol).  The  results  also  indicate  that  the  transition 
state  for  reaction  (I)  has  three-fold  symmetry.  Evidence  of  a  step-wise  decomposition  reaction, 

sym-triazine  H2C2N2  +  HCN  ->  3  HCN,  (II) 

was  not  observed.  In  the  248-nm  experiments,  the  energy  measured  in  the  product  translational  and 
internal  modes  is  consistent  with  an  equipardtioning  of  the  available  product  energy.  The  193-nm 
experiments,  however,  showed  a  nonstatistical  product  energy  distribution  with  99%  of  the  available 
product  energy  partitioned  into  the  internal  modes  of  the  HCN  products  (Ondrey  and  Bersohn  1984).  A 
time-resolved  infrared  fluorescence  study  of  the  HCN  formed  from  193-nm  jrfiotolysis  of  sym-triazine 
showed  that  the  HCN  bending  vibrations  are  excited  preferentially,  with  "no  evidence  of  a  hot  rotational 
population"  (Goates,  Chu,  and  Flynn  1984).  A  simple  harmonic  oscillator  analysis  of  the  products 
indicates  that  the  bending  quanta  of  HCN  formed  from  photolysis  of  sym-triazine  at  193  nm  is  70  times 
larger  than  the  number  of  C-H  stretching  quanta  (Goates,  Chu,  and  Flynn  1984). 

The  anmnalous  difference  in  the  product  energy  distributions  for  the  248-  and  193-nm 
photodissociation  experiments  was  attributed  to  transitions  to  different  excited  states  that  subsequently 
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cross  to  the  ground  state  potential  energy  surface  at  different  regions  of  configuration  space  (Migidichian 
1947).  The  geometric  structures,  levels,  and  types  of  vibrational  excitation  in  the  ground-state  molecule 
are  dependent  on  the  region  of  configuration  space  at  which  the  crossings  occur.  If  the  transition  results 
in  a  vibrationally  excited  and  highly  distorted  ground-state  molecule  that  reacts  before  vibrational  or 
structural  relaxation  occurs,  a  nonstatistical  partitioning  of  energy  could  result  in  the  product  molecules. 
This  could  explain  the  differences  in  product  energy  distributimis  for  the  two  i*otolysis  energies  (Ondrey 
and  Bersohn  1984). 

The  experiments  clearly  show  that  decomposition  occurs  in  a  concerted  marmer  under  193-  and 
248-nm  i^tolysis  (Ondrey  and  Bersohn  1984;  Goates,  Chu,  and  Flytm  1984).  Concerted  triple 
dissociation  reactions  are  uncommon,  but  not  unknown.  There  are  at  least  two  other  cyclic  molecular 
systems  that  have  also  been  shown  to  decompose  through  a  concerted  triple  dissociation  mechanism, 
5-tetrazine  (King  et  al  1977;  Zhao  et  al.  1989)  and  2,4,6-hexahydro-l,3,5-trinitro-l,3,5-triazine  (Zhao, 
Hintsa,  and  Lee  1988). 

Concerted  triple  association  of  three  HCN  molecules  to  form  sywj-triazine  seems  more  improbable  than 
the  reverse  reaction,  the  concerted  triple  dissociation  in  reaction  (I).  In  a  triple  association  reaction, 
entropic  effects  associated  with  bringing  three  HCN  molecules  together  in  a  cmicerted  fashion,  especially 
in  the  gas  i^rase,  would  seem  so  large  as  to  prohibit  reaction.  However,  entropic  hindrance  to  association 
would  be  reduced  if  a  prereaction  intermediate  with  a  structure  favorable  to  concerted  triple  associ^on 
could  be  easily  formed.  Experiment  shows  that  such  a  species  exists  (Jucks  and  Miller  1988).  It  is  a 
cyclic  hydrogen-bonded  cluster  (HCN)3.  RotationaUy  resolved  spectra  of  the  infrmed  active  doubly 
degenerate  C-H  stretch  for  the  trimer  cluster  is  consistent  with  an  oblate  planar  symmetric  top  (Jucks  arid 
Miller  1988).  Unfortunately,  the  experimental  results  could  not  provide  details  of  this  cyclic  structure. 
Subsequent  ab  initio  calculations  on  hydrogen-bonded  (HCN)3  clusters  confirmed  that  the  cyclic  cluster 
exists  in  conjunction  with  linear  HCN  trimer  chains  (Kumig,  Lischka,  and  Karpfen  1990).  Zero-point 
corrected  relative  ertergies  of  the  two  clusters  differ  by  only  0.5  kcalAnol  (Kurnig,  Lischka,  and  Karpfen 
1990).  It  is  conceivable  that  HCN  gas,  under  containment,  would  form  clusters,  including  the 
hydrogen-bonded  cyclic  trimer.  The  HCN  molecules,  weakly  bound  in  such  a  cyclic  cluster,  would  then 
be  in  a  sterically  favorable  arrangement  for  a  concerted  association: 

3  HCN  (HCN)3  -¥  sym-triazine.  (Ill) 

Reaction  energy  ai^ropriately  imparted  to  the  cluster  would  allow  reaction  (III)  to  occur. 
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The  focus  of  this  study  is  to  investigate  reaction  mechanisms  and  eneigetics  for  the  decomposition  and 
formation  of  sym-triazine  [reactions  (I)-(ni)].  We  also  investigate  vibrational  coupling  of  ^m-triazine 
and  the  (HCN)3  cluster  with  the  reaction  coordinate  for  reactions  (I)  and  (III).  Our  results  also  examine 
basis  set  dependence  for  die  calculations,  as  weU  as  effects  of  increased  electron  correlation. 

2.  METHODS 

All  calculations  reported  in  this  woilc  were  done  using  the  Gaussian  94  set  of  programs  (Frisch  et  al. 
1995).  Critical  points  on  the  sym-triazine  potential  energy  surface  corresponding  to  reactions  (I)-(ni)  were 
determined  through  MP2  geometry  optimizations  and  characterized  through  normal  mode  analyses  using 
the  6-31 G**  (Hehre,  Ditchfield.  and  Pople  1972;  Hariharan  and  Pople  1973;  Gordon  1980),  6-311-h-G** 
(McLean  and  Chandler  1980;  Krishnan  et  al.  1980),  and  cc-pVTZ  (Woon  and  Dunning  1993;  Kendall, 
Dunning,  and  Harrison  1972;  Dunning  1989)  basis  sets.  QCISD(T)  energy  refinements  were  calculated 
for  each  MP2-optimized  structure  with  its  corresponding  basis  set.  All  calculations  were  run  with  frozen 
core  orbitals,  and  die  geometry  optimizations  met  the  default  convergence  criteria  given  by  Frisch  et  al. 
(1995).  Intrinsic  reaction  coordinate  (IRC)  calculations  were  performed  at  the  MP2/6-31G**  level.  The 
IRC  calculation  leading  from  the  transition  state  structure  of  reaction  (I)  towards  isolated  HCN  confirmed 
that  this  transition  state  is  also  the  saddle  point  for  reaction  (III).  All  IRC  calculations  proceeded  to 
minima  as  determined  by  the  default  convergence  criteria  in  Frisch  et  al.  (1995).  Saddle  points  for 
reaction  (II)  were  not  determined  for  reasons  discussed  as  follows. 

3.  RESULTS  AND  DISCUSSION 

Molecular  stmctures  for  jym-triazine,  the  dimer  H2C2N2,  the  cyclic  (HCN)3  cluster,  the  HCN 
molecule,  and  the  transition  state  for  reactions  (I)  and  (III)  are  shown  in  Figure  1.  Subsequent  normal 
mode  analyses  characterized  each  species.  The  molecular  properties  are  shown  in  Table  1.  The  atom 
labels  in  Table  1  are  consistent  with  the  structures  shown  in  Figure  1,  and  the  internal  coordinates  used 
in  our  discussion  are  shown  in  the  transition  state  structure  in  Figure  1(b).  Although  the  atoms  are  not 
labeled  on  the  transition  state  structure  in  Figure  1(b),  the  labeling  of  the  atoms  follows  the  same  pattern 
around  the  ring  as  in  Figure  1(a).  The  internal  coordinates  shown  in  Figure  1(b)  will  be  discussed  next 
and  are  the  same  for  all  structures  illustrated  in  Figure  1. 
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c)  the  stable  dimer  species  H 


with  the  stepwise  dissociation  miechanism  Ireaction  (IDI:  (d)  HCN:  and  (e~)  the  hydrogen-bonded 


(HCN)^  cluster  located  on  the  concerted  triple  association/dissociation  path. 
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•  Energies  in  kcalAnol,  cJisUnces  in  A.  bond  angles  in  degrees,  and  frequencies  in  cin‘^  ^  Lancaster  and  Stoicheff  1956.  ®  Lancaster,  Stamm,  and  Cdthup  1961.  ^  Hcrrf)eig  1945;  Huber  and 
Herzberg  1979.  ®  Kumig,  lischka,  and  Karpfen  1990.  ^  Atom  numbering  illustrated  in  Figure  1.  *  Experimental  value  is  3273.545  cm“'  (Jucks  and  Miller  1988). 


3.1  Geometries.  All  calculated  structural  parameters  of  s)wi-triazine  are  within  1%  of  the 
experimental  values,  lowing  very  little  variation  with  the  size  of  the  basis  set.  The  HCN  geometries  are 
in  good  agreement  with  experiment  for  all  the  basis  sets,  with  the  errors  ranging  from  <2%  (6-310**)  to 
<1%  (cc-pVTZ).  The  structure  of  the  (HCN)3  cluster  has  not  been  determined  experimentally,  although 
Jucks  and  Miller  (1988)  suggest  two  structures  that  are  consistent  with  the  rotational  constants  determined 
from  the  rovibrational  spectnim  of  the  C-H  asymmetric  stretch.  The  intermolecular  C-C  distances  for  the 
two  structures  suggested  by  Jucks  and  Miller  are  3.62  and  3.83  A,  respectively.  The  C-C  intemuclear 
distances  of  the  (HCN)3  clusters  calculated  using  the  6-31G**,  6-31 1++G**,  and  cc-pVTZ  basis  sets  are 
3.63, 3.66,  and  3.60  A,  respectively.  These  distances  are  consistent  with  those  of  Structure  A  in  Figure  5 
of  Jucks  and  Miller  (1988).  Ab  initio  calculations  using  the  average  coupled  pair  approximation  (ACPF) 
method  and  extended  basis  sets  (Kumig,  Lischka,  and  Karpfen  1990)  predicted  the  cyclic  trimer  with  N-H 
intermolecular  distances  of  2.5532  A.  Comparison  of  our  results  with  the  ACPF  calculations  is  shown 
in  Table  1.  The  structures  from  the  two  methods  are  similar.  The  cyclic  cluster,  illustrated  in  Figure  1(e), 
consists  of  three  linear  HCN  molecules  with  H-N  distances  of  2.5  A. 

Structures  and  energies  along  the  reaction  coordinate  determined  from  MP2/6-31G**  IRC  calculaticms 
are  shown  in  Figure  2.  The  structures  along  the  reaction  path  have  three-fold  symmetry.  Because  of  this, 
there  are  six  unique  internal  coordinates  that  describe  the  geometric  changes  along  the  reaction  path  and 
are  defined  in  Figure  1(b).  These  unique  internal  coordinates  include  two  types  of  C-N  bonds:  one  that 
becomes  intermolecular  in  the  (HChOs  cluster  (denoted  CN2),  and  one  that  remains  an  intramolecular  bond 
along  the  reaction  path  (denoted  (3^1).  The  C-H  distances  and  HCN,  NCN,  and  CNC  angles  make  up 
the  rest  of  the  six  internal  coordinates.  Changes  of  these  internal  coordinates  along  the  reaction  pafli  are 
shown  in  Figure  3.  Negative  values  of  the  reaction  coordinate  correspond  to  the  (HCN)3  cluster  region 
of  the  PES,  and  posative  values  along  the  reaction  coordinate  correspond  to  the  sym-triazine  region  of  the 
PES.  The  reaction  path  coordinate  value  0.0  corresponds  to  the  transition  state  connecting  the  (HCN)3 
cluster  and  ^/n-triazine  minima.  The  internal  coordinates  that  change  die  most  along  this  reaction  path 
are  the  CN2  bond  length  and  the  HCN  angle.  The  CN2  bond  increases  almost  linearly  as  the  jym-triazine 
dissociation  progresses  to  the  (HCN)3  cluster.  The  other  CN  bond  type,  CNl,  becomes  shorter  as 
sym-triazine  decomposes,  and  almost  equals  the  equilibrium  CN  bond  distance  in  the  (HCN)3  cluster  and 
isolated  HCN  by  the  time  it  reaches  the  saddle  point.  The  CNl  bond  distance  changes  only  slightly  as 
the  (HCN)3  minimum  is  approached  after  crossing  the  saddle  point.  The  C-H  bond  distances  do  not 
change  significantly  during  the  reaction.  The  HCN  angle  increases  almost  monotonically  as  tire  (HCN)3 
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Figure  2.  MP2/6-31G**  reaction  oath  from  the  IRC  calculation  for  reactions  (I)  and  (III).  Every  tenth 
point  of  the  IRC  is  shown.  In  addition  to  the  stable  and  transition  state  structures,  three  other 
structures  alone  the  reaction  oatti  have  been  shown  to  enable  the  reader  to  visualize  the 
mechanism  of  concerted  triple  association  and  dissociatioa  The  (HCN)^  is  illustrated  in  the 
far-left  structure  of  the  fieure.  The  5v/w-triazine  molecule  is  represented  bv  the  far-right 


CNl  (A) 


(amu^^^-bohr)  (amu^'^^-bohr) 


Figure  3.  Internal  coordinate  changes  alone  the  MP2/6-31G**  IRC  for  the  conceited  triple  association/ 
dissociation  reactions.  For  definition  of  internal  coordinates,  see  Figure  1(c). 


8 


HCN  (degrees)  NCN  (degrees)  CNC  (degrees) 


cluster  is  ^roached.  This  angle  varies  from  117°  in  ^ym-triazine  to  139°  at  the  transition  state.  It 
continues  to  open  up  to  180°  in  the  (HCN)3  cluster  region.  The  ring  angles,  NCN  and  CNC,  do  not  show 
monotonic  behavior.  The  NCN  angle  first  decreases  alter  crossing  the  transition  state  and  then  opens  up 
rt^idly  as  the  molecule  ^proaches  the  cluster  minimum.  The  changes  in  the  CNC  angle  along  the 
reaction  path  behave  in  an  almost  equal  but  opposite  manner  to  the  changes  in  the  NCN  angles. 

Attempts  were  made  to  locate  transition  states  leading  to  dimer  formation  and  decomposition  and 
stable  structures  starting  from  various  dimer  configurations.  The  only  structure  located  corresponds  to  a 
local  minimum  on  the  PES.  Its  structure  is  that  of  a  parallelogram  of  die  heavy  atoms  with  internal  angles 
for  CNC  and  NCN  of  100®  and  79°,  respectively  [Figure  1(c)].  The  HCN  angles  are  far  from  linear  in 
this  structure  (see  Table  1). 

3.2  Frequencies.  MP2  harmonic  vibrational  frequencies  of  the  critical  points  described  previously 
were  calculated  using  the  aforementioned  basis  sets  and  are  shown  in  Table  1.  MP2/cc-pVTZ  eigenvectors 
and  corresponding  harmonic  fiequencies  of  all  critical  points  (except  HCN)  for  reactions  (I)  and  (III)  are 
shown  in  Figures  4-6.  The  predicted  fiequencies  of  HCN  agree  with  experiment  to  within  7%  for  all 
basis  sets.  Predictions  of  sym-triazine  frequencies  are  within  7.5%  of  experimental  values  for  aU  basis 
sets.  There  is  only  one  vibration  that  has  been  measured  and  assigned  for  the  (HCN)3  cluster.  It 
corresponds  to  the  doubly  degenerate  C-H  asymmetric  stretch  at  3274  cm~^  and  contains  rotationally 
resolved  bands  (Jucks  and  Miller  1988).  Our  value  at  3438  cm~^  is  5%  higher  than  the  experimental 
value.  Judes  and  MiUer  provide  rotational  constants  A"=B"=2C"=0.0822  cm~^  deteimined  from  fits  to 
the  rotaticmal  transitions  for  a  planar  oblate  symmetric  top  (Jucks  and  Miller  1988).  Our  calculated  values 
using  the  6-31G**,  6-31 1++G**,  and  cc-pVTZ  basis  sets  are  0.0827,  0.0821,  and  0.0840  cm~^ 
respectively. 

Lancaster,  Stamm,  and  Colthup  (1961)  reported  the  infrared  and  Raman  spectra  for  sy/n-triazine  and 
made  band  assignments  based  upon  a  normal  coordinate  analysis.  Using  the  description  of  the  modes 
given  by  Lancaster,  Stamm,  and  Colthup,  we  were  able  to  match  all  vibrational  modes  to  experimental 
assignments  except  for  the  experimental  modes  V4  (1617  cm“^),  Vg  (1410  cm“^),  and  V9  (1174  cm“^). 
Lancaster,  Stamm,  and  Colthup  predicted  doubly  degenerate  modes  Vg  and  V9  should  correspond  to  ring 
and  C-H  rocking  vibrations,  respectively.  If  we  assume  that  the  MP2/cc-pVTZ  modes  with  frequency 
1449  cm“^  can  be  assigned  to  experimental  Vg  and  the  modes  with  frequency  1 199  cm“^  to  experimental 


9 


Each  eigenvector  is  labeled  with  its  cone! 


inding  fteouencv  (in  cm 


Vg,  then  the  ab  initio  predictions  agree  with  experiment  to  within  3%.  Visual  inspection  of  the 
MP2/cc-pVTZ  eigenvectors  corresponding  to  ab  initio  frequencies  1449  and  1 199  cm“^  indicate  that  these 
two  sets  of  degenerate  vibrations  cannot  be  distinctly  classified  as  C-H  rocks  or  ring  vibrations.  The 
remaining  ab  initio  eigenvector  that  we  could  not  assign  to  anything  else  (1405  cm~^)  was  assigned  to 
experimental  mode  V4  with  13%  disagreement  However,  mode  V4  is  one  of  the  two  A  2  fundamentals 
that  are  inactive  and  were  not  observed  (lucks  and  Miller  1988).  The  frequency  for  this  fundamental  that 
was  reported  by  Lancaster,  Stamm,  and  Colthup  (1961)  was  calculated  from  estimated  force  constants  and 
considered  "very  approximate." 

The  zero-point  eneipes  for  sym-triazine  with  die  6-31G**,  6-31 1++G**,  and  the  cc-pVTZ  basis  sets 
are  41.6, 40.9,  and  41.2  kcal/mol,  respectively,  and  compare  well  with  experimental  zero-point  energy  of 
40.5  kcal/mol  (Lancaster,  Stamm,  and  Colthup  1961).  The  HCN  zero-point  energies  for  the  three  basis 
sets  (ranging  from  9.9  to  10.0  kcal/mol)  agree  weU  with  the  experimental  value  (9.8  kcalAnol)  (Herzberg 
1945;  Huber  and  Hereberg  1979). 

3.3  Vibrational  Coupling.  It  has  been  shown  that  projection  of  the  eigenvector  of  a  vibrational  mode 
onto  the  eigenvector  corresponding  to  the  direction  along  the  reaction  path  is  related  to  the  coupling  of 
that  vibrational  mode  with  the  reaction  coordinate  (Waite  and  Miller  1981;  Rice,  Grosh,  and  Thomp)Son 
1995).  Waite  and  Miller  showed  that  the  unimolecular  decay  rates  of  the  Henon-Heiles  model  behave 
statistically  for  all  energies,  even  though  this  system  has  quasiporiodic  classical  motion  at  low  eneigies, 
conditions  under  which  mode  sp)ecificity  might  be  exp)ected  (Waite  and  Miller  1981).  They  suggested  that 
the  statistical  behavior  was  due  to  coupling  of  the  intramolecular  motions  of  the  model  with  the 
dissociative  reaction  coordinates  for  the  system,  since  there  was  some  degree  of  projection  of  all  of  the 
vibrational  modes  onto  the  reaction  coordinates  for  the  Henon-Heiles  model.  In  other  words,  there  were 
no  vibrational  modes  in  which  energy  could  be  trapp}ed.  To  further  investigate  the  role  of  vibrational 
coupling  with  the  reaction  coordinate,  they  modified  the  Henon-Heiles  px>tential  such  that  vibrational 
modes  exist  that  do  not  project  onto  the  reaction  coordinate  (Waite  and  Miller  1981).  Mode  specificity 
in  the  decmnpxisition  was  then  induced.  Rice,  Grosh,  and  Thompson  (1995)  showed  that  rates  of 
unimolecular  decompx>sition  through  competing  channels  were  enhanced  and  branching  ratios  changed 
upxin  excitation  of  vibrational  modes  of  a  reactant  that  project  strongly  onto  the  reaction  coordinates  of 
the  system  (Rice,  Gro^,  and  Thompson  1995).  Our  goal  in  this  section  is  not  to  investigate  mode-spieciric 
dynamics  for  this  system.  Rather  we  wish  to  determine  which  of  the  intramolecular  motions  of 
jym-triazine  and  the  (HCN)3  cluster  couple  most  strongly  with  the  reaction  coordinate  for  reactions  (I)  and 
(III).  This  will  suggest  mechanisms  for  enei^gy  transfer  from  excited  vibrational  modes  to  the  reaction 
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coordinate.  The  studies  by  Waite  and  MiUer  (1981)  and  Rice,  Grosh,  and  Thompson  (1995)  have  clearly 
shown  a  correlation  between  projectirai  of  vibrational  modes  onto  tiie  reaction  coordinate  and  coupling 
with  the  reactiwi  coordinate.  Thus,  similar  analyses  based  on  projections  of  the  vibrational  modes  of 
jyw-tiiazine  and  the  (HCNOj  cluster  onto  the  reaction  coordinate  should  verify  whether  the  intramolecular 
motions  of  diese  molecules  are  coupled  with  the  reaction  path. 

As  set  forth  in  Rice,  Grosh,  and  Thompson  (1995),  we  calculated  local  normal  modes  (Miller,  Handy, 
and  Adams  1980)  for  points  along  the  reaction  path  for  reaction  (HI)  at  the  MP2/6-31G**  level.  The 
inhniteamal  translations  and  rotations  were  projected  out,  leaving  3N-7  vibrational  modes  of  the  molecule 
and  the  eigenvector  corresponding  to  the  direction  along  the  reaction  path  (Miller,  Handy,  and  Adams 
1980),  We  projected  the  eigenvectors  corresponding  to  the  harmonic  vibrational  frequencies  of 
equilibrium  jy/n-triazine  and  the  (HCN)3  cluster  onto  the  eigenvector  associated  with  the  direction  along 
the  reaction  path  (Miller,  Handy,  and  Adams  1980)  for  selected  reaction  coordinate  values.  The  results 
of  the  projections  of  5yw-triazine  and  the  (HCN)3  cluster  vibrational  eigenvectors  onto  the  reaction  path 
eigenvectors  are  shown  in  Rgures  7  and  8,  respectively.  Only  those  modes  that  have  projections  greater 
than  0.05  are  shown  in  these  figures.  There  are  three  vibrational  modes  of  aym-triazine  that  project 
strongly  all  along  the  reaction  path  fpr  reactions  (I)  and  (III);  they  correspond  to  1015,  1149,  and  1430 
cm“*.  These  modes  correspond  to  two  ring  symmetric  breathing  modes  and  an  HCN  symmetric  bending 
mode,  respectively.  Additionally,  a  fourth  vibrational  mode  (1256  cm  ^)  jnnjects  onto  portions  of  the 
reaction  path,  but  its  projection  is  not  as  strong  as  the  other  vibrational  modes.  Its  motion  can  best  be 
described  as  being  the  most  similar  to  the  eigenvector  associated  with  the  imaginary  frequency  at  the 
saddle  point  (see  Figure  5).  There  are  two  vibrational  modes  for  the  (HCT^j  cluster  that  project  strongly 
onto  the  reaction  path;  these  are  both  symmetric  ring-breathing  modes  (105  and  149  cm'^).  Two  other 
modes  project  less  strwigly  onto  the  reaction  path  (740  and  3496  cm“^);  both  consist  mainly  of  hydrogenic 
motions. 

It  is  clear  that  there  are  intramolecular  motions  of  ay/n-triazine  and  the  (HCN)3  cluster  that  couple 
strongly  with  the  reaction  coordinate,  indicating  the  most  likely  routes  through  which  reaction  energy 
distributed  in  either  of  these  species  can  transfer  efficiently  to  the  reaction  coordinate  for  reactions  (I)  and 
(HI). 

3.4  Energetics.  The  relative  and  absolute  energies  of  each  critical  point  calculated  at  the  MP2  and 
QCISDCO  levels  with  the  three  different  basis  sets  for  ^m-triazine,  H2C2N2,  the  (HCN)3  cluster,  HCN, 
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Figure  8.  Projection  of  the  eigenvectors  corresponding  to  the  MP2/6-31G**  haimonic  vibrational  freQuendes  of  the  (HCK)^  cluster  onto 
eigenvectors  for  points  on  the  reaction  path  that  correspond  to  the  direction  along  the  oath  for  the  concerted  trinle  as!«iriatinf^ 
decomposition  reactions  freactions  (I)  and  (IIDI.  Eigenvectors  with  projections  greater  than  0.05  are  illustrated,  and  labeled  with  the 
corresponding  MP2/6-31G**  harmonic  vibrational  frequency  (in  cm"^)' 


and  the  saddlepoint  are  listed  in  Table  2.  (A  comparison  of  the  zero-point-energy-corrected  MP2  barriers 
to  fonnation  of  sym-triazine  from  the  (HO^)3  cluster  and  to  decomposition  of  iym-triazine  with  QC1SD(T) 
for  all  the  basis  sets  is  shown  in  Figure  9.)  The  MP2/6-31G**  and  MP2/6-31 1++G**  calculations  predict 
similar  barrier  heights  for  the  formation  of  iym-triazine  from  the  (HCN)3  cluster  (62.2  and  62.4  kcalAnol, 
respectively).  The  MP2/cc-pVTZ  prediction  of  this  barrier  is  lower  (57.4  kcalAnol).  QQSEKD  energy 
refinements  for  the  6-31G**  and  6-31 1++G**  barriers  decrease  the  MP2  results  slightly  (61.1  and  61.7 
kcal/mol,  respectively).  The  QCISD(T)  refinement  of  the  MP2/cc-pVTZ  barrier,  however,  increases  the 
MP2  result  by  0.7  kcalAnol  (58.1  kcal/mol).  Thermal  activation  barriers  for  this  reaction  have  not  been 
determined  experimentaUy;  thus,  we  caimot  gauge  the  accuracy  of  these  barrier  heights.  The  only 
information  about  the  barrier  height  for  the  formation  of  jym-triazine  from  isolated  HCN  is  an  upper  limit 
of  72  kcalAnol  (Ondrey  and  Bersohn  1984).  Our  values  are  well  under  this  limit. 

The  zero-point-corrected  MP2  barriers  to  decomposition  of  aym-triazine  for  the  6-31G**,  6-31 1++G**, 
and  cc-pVTZ  basis  sets  are  80.7,  77.7,  and  79.6  kcal/mol,  respectively.  QQSD(T)  refinements  of  the 
barriers  for  these  three  basis  sets  are  81.6,  79.0,  and  81.2  kcalAnol,  respectively;  well  below  the 
experimental  upper  limit  of  1 15  kcal/mol  (Ondrey  and  Bersohn  1984).  QaSD(T)  energy  refinements  did 
not  significantly  change  the  MP2  results.  Both  reactions,  however,  showed  some  sensitivity  to  basis  set 
size. 

Experimental  and  calculated  reaction  enthalpies  (corrected  to  T-298  K)  for  reaction  (I)  are  given  in 
Table  3.  The  calculated  reaction  enthalpies  for  reaction  (I)  range  from  27.8  to  35.5  kcal/mol,  at  least  18% 
lower  than  the  experimental  value  of  43.2  kcalAnol  (Ondrey  and  Bersohn  1984).  The  QQSD(T)//MP2/ 
cc-pVTZ  value  of  35.5  kcalAnol  is  in  closest  agreement  with  experiment.  The  fonnation  enthalpies  of  the 
(HCN)3  cluster  and  H2C2^2  ^  ^  shown  in  Table  3. 

The  zero-point-corrected  MP2  energies  of  the  H2C2N2  [Figure  1(c)]  plus  HCN  relative  to  isolated 
HCN  are  63.8,  64.0,  and  61.3  kcalAnol  for  the  6-31G**,  6311-h-G**,  and  cc-pVTZ  basis  sets, 
respectively.  Corresponding  zero-point-corrected  QC1SD(T)  values  are  57.4,  57.6,  and  55.9  kcalAnol. 
These  values  are  higher  than  the  zero-point-eneigy-corrected  barriers  to  concerted  triple  association  at  all 
levels  for  each  basis  set,  with  a  difference  of  6.7  kcal/mol  at  the  highest  level  of  theory.  Figure  10 
illustrates  this  for  the  QClSD(T)/cc-pVTZ  results.  Even  if  we  assume  that  there  are  no  barriers  to 
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Table  2.  Absolute  Energies  (Hartrees)  and  Relative  Energies  With  and  Without  JZero-Point  Energy  (ZPE)  Corrections  (kcal/mol)  of 
Species  on  the  sym-Triazine  Potential  Energy  Surface 
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Table  3.  Temperature-Corrected  (T-298  K)  Enthalpies  (kcal/mol) 


Level  of  Theory 

Reaction  Endotheimicity 

sym-triazine  3HCN 

Experiment* 

43.2 

MP2//MP2/6-31G** 

31.5 

(K:ISD(T)//MP2/6-3  IG** 

33.0 

MP2//MP2/6-311++G** 

27.8 

QaSD(T)//MP2/631 1++G** 

29.5 

MP2//MP2/CC-PVTZ 

35.2 

(y:iSD(T)//MP2/cc-pVTZ 

35.5 

3HCN  ^  (HCN)3 

MP2//MP2/6-31G** 

-9.4 

QaSD(D//MP2/6-31G** 

-8.9 

MP2//MP2/6-311++G** 

-8.7 

(JCISD(T)//MP2/631  1-h-G** 

-8.4 

MP2//MP2/cc-pVTZ 

-9.3 

QCISD(T)//MP2/cc-pVTZ 

-8.7 

3HCN  ->  H2C2N2  +  HCN 

MP2//MP2/6-31G** 

62.1 

QCISD(T)//MP2/6-31G** 

55.7 

MP2//MP2/6-311++G** 

62.5 

QCISD(T)//MP2/631 1++G** 

56.1 

MP2//MP2/cc-pVTZ 

59.7 

QCISD(T)//MP2/cc-pVTZ 

54.3 

^  Ondrey  and  Bersohn  1984. 


formation  of  H2C2N2  or  the  rrecessary  HCN  insertion  into  the  H2C2^2  molecule  to  form  iym-triazine 
[reverse  of  reaction  (II)],  comparing  the  relative  energies  of  the  H2C2N2  +  HCN  minimum  and  the  saddle 
point  for  reaction  (HI)  indicates  that  the  reaction  (II)  pathway  is  higher  in  energy  than  that  of  the  conceited 
triple  association  reaction  (III).  Because  we  are  interested  mainly  in  the  low-energy  pathway  for  the 
formation/decomposition  reactions  of  sy/n-triazine,  we  did  not  investigate  this  reaction  path  further.  We 
attempted  to  locate  additional  dimer  structures,  but  none  were  found. 
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Figure  10.  Schematic  of  the  overall  energetics  for  the  reactions  3HCN  ->  jv/w-triazine.  Zero-ooint- 
conected  energies  at  the  OCISD(D//MP2/cc-dVTZ  level  are  listed  in  kcal/mole.  The  higher 
energy  pathway  is  the  stepwise  association  mechanism  freverse  of  reaction  (11)1:  the  lower 
energy  path  is  the  concerted  triple  association  mechanism  fieaction  (PI. 
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These  results  support  the  experimental  observations  for  reaction  (I)  (Ondrey  and  Bersohn  1984; 
Goates,  Chu,  and  Flynn  1984).  First,  the  zero-point-energy-corrected  barrier  for  .^m-triazine 
decomposition  [reaction  (I)]  is  within  the  experimentally  determined  upper  limit.  Secondly,  the 
low-energy  decomposition  reaction  is  a  concerted  triple  dissociation  rather  than  a  step-wise  decomposition 
reaction  sudi  as  reaction  (II),  in  agreement  with  the  i*otodissociation  ejqwrimental  results.  Additionally, 
the  structures  of  species  along  the  reaction  path  for  this  process  could  explain  the  observed  nonstatistical 
energy  distribution  in  the  bending  mode  over  the  stretches  (Goates,  Chu,  and  Flyrm  1984).  The  saddle 
point  structure  for  this  reaction  is  similar  to  that  of  sym-triazine.  The  HQ^  angles  in  this  structure  are 
140^.  Thus,  (mce  this  barrier  is  overcome  and  the  system  is  in  the  HCN  region  of  the  potential  energy 
surface,  the  HCN  angles  must  opai  up  40“  to  reach  both  the  (HChOj  cluster  and  HCN  product  values  of 
180“.  The  large  rhangt*  in  the  HCN  angle  of  the  transition  state  species  to  the  clusterdsolated  HCN 
molecules  could  be  the  source  of  the  large  HCN  bending  excitation  that  was  observed  by  Goates,  Chu, 
and  Flynn  (1984)  in  their  study  of  the  jwoduct  energy  distributions  of  sym-triazine  photodissociation  at 
193  run.  The  CN  and  CH  bond  distances,  on  the  other  hand,  differ  at  most  by  5%  between  the  saddle 
point  and  cyclic  (HCN)3  cluster  geometries.  Thus,  the  most  significant  geometric  changes  in  these  internal 
coordinates  have  occurred  before  the  saddle  point  is  reached  and  could  explain  why  vibrations 
corresponding  to  CN  and  CH  stretches  were  not  excited  upon  decomposition. 

Our  reaction  path  calculations  indipate  that  the  weakly  bound  cyclic  (HCN)3  cluster  is  a  reaction 
intermediate  on  the  sym-triazine  decomposition  pathway.  However,  the  energy  available  to  products  upon 
crossing  the  reaction  barrier  is  so  large  that  this  reaction  intermediate  would  not  be  long-lived  and 
probably  not  detectable  during  the  decomposition  process.  Our  results  clearly  show  that  the  low-energy 
path  to  dissociation  of  jym-triazine  is  throu^  a  concerted  triple  dissociation.  We  were  unable  to  locate 
a  saddle  point  for  decomposition  of  the  (HCN)3  cluster  to  form  isolated  HCN.  Thus,  we  are  assuming  that 
the  barrier  to  (HCN)3  decomposition  toward  isolated  HCN  is  merely  the  endothermicity  of  the  cluster 
relative  to  isolated  HCN. 

4.  CONCLUSIONS 

We  have  presented  an  ab  initio  study  of  formation  and  decomposition  reactions  of  iy/n-triazine.  Two 
decomposition  pathways  were  examined:  a  concerted  triple  decomposition  reaction 

sy/n-triazine  — >  3  HCN  (I) 
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and  a  st^wise  decomposition  reaction 


s}'m-triazine  ->  H2C2N2  ^  HOi,  (II) 

Our  best  estimate  of  the  eneigy  required  for  reactions  (I)  and  (11)  are  81.2  and  87.0  kcai/mol  [QCISDCD// 
MP2/cc*pVTZ],  respectively.  These  predicted  enei:gy  requirements  for  reactions  (I)  and  (II)  are  well  below 
the  experimentally  determined  upper  bound  (115  kcal/mol)^. 

IRC  calculations  starting  from  the  transition  state  of  reaction  (I)  toward  isolated  HCN  led  to  a  local 
minimum  on  the  potential  eneigy  surface  that  corresponds  to  a  weakly  bound  cyclic  (HCN)3  cluster. 
Thus,  the  reaction  path  connects  the  sym-triazine  minimum  with  the  cyclic  (HCN)3  cluster,  a  reaction 
intermediate  leading  to  isolated  HCN.  Our  best  estimate  of  its  eneigy  relative  to  isolated  HCK  is  -8.9 
kcal/mol  [QCISD(T)//MP2/cc-pVTZ].  The  reverse  of  reaction  (I),  the  association  of  HCN  molecules  to 
fonn  sym-triazine,  is  therefore 


3  HCN  -¥  (HCN)3  -» jym-triazine.  (HI) 

It  is  likely  that  fonnation  of  the  (HCN)3  cluster  is  a  critical  step  in  the  association  reaction  (III)  as  it 
removes  significant  steric  hindrance  to  the  concerted  triple  association  of  the  HCN  molecules  to  form 
jym-triazine.  However,  this  cluster  is  probably  very  short  lived  upon  decomposition  of  sym-triazine,  due 
to  the  excess  eneigy  available  to  products  upon  crossing  the  saddle  point  for  reaction  (I). 

The  zero-point-eneigy-corrected  reaction  barrier  to  triple  association  [reaction  (III)]  is  lower  than  the 
eneigy  needed  to  form  the  HCN  +  H2C2N2,  the  stable  intermediate  complex  for  the  step-wise  reaction  for 
formation  or  decomposition  of  sym-triazine  for  all  basis  sets  and  all  levels  of  electron  correlation.  Our 
best  estimate  of  the  eneigy  difference  between  HCN  -f  H2C2N2  and  the  barrier  to  triple  association 
(relative  to  isolated  HCN)  is  6.7  kcal/mol  [QCISD(T)//MP2/cc-pVTZ].  Any  energy  barrier  to  formation 
of  H2C2N2  or  insertion  of  a  third  HCN  into  this  intermediate  to  form  jy/n-triazine  cannot  be  lower  than 
the  oiergy  of  this  intermediate.  Thus,  a  step-wise  reaction  for  formation  or  decomposition  of  jym-triazine 
is  eliminated  as  the  low-eneigy  path  in  favor  of  the  conceited  route. 

Critical  points  m  the  sym-triazine  potential  energy  surface  were  located  through  MP2  geometry 
optimizations  using  the  6-31G**,  6-31 1++G**,  and  cc-pVTZ  basis  sets,  and  characterized  through  normal 
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mode  analyses.  Energy  refinements  for  each  point  were  done  at  the  QCISD(T)  level  using  the  same  three 
basis  sets.  It  was  shown  diat  relative  energy  changes  between  the  two  levels  of  electron  correlation  were 
small.  Structural  changes  for  the  critical  points  with  increasing  basis  set  were  also  small,  and  all  predicted 
structures  were  in  good  agreement  with  experiment,  where  available.  Because  the  geometries  of  HCN  arul 
sym-triazine  are  described  extremely  well  at  the  MP2  level  with  all  three  basis  sets,  this  suggests  accurate 
geometries  are  available  without  requiring  large  basis  sets.  The  MP2  frequencies  show  reasonable 
comparison  with  experimental  values  where  available.  Temperature-corrected  (T=298  K)  reaction 
enthalpies  for  dissociation  of  syw-triazine  [reaction  (I)],  however,  were  lower  than  experiment  (Ondrey 
and  Bersohn  1984)  by  18%. 

IRC  calculations  leading  frcxn  the  transition  state  for  formation/dissociation  of  sym-triazine  provided 
mechanistic  insight  into  these  reactions.  The  species  along  the  reaction  path  have  three-fold  symmetry, 
in  agreement  with  photodissociation  experiments  that  clearly  show  concerted  triple  dissociation  of  the 
syw-triazine  (Ondrey  and  Bersohn  1984).  The  large  change  in  the  HCN  angle  going  from  the  transition 
state  to  the  (HCN)3  cluster/HCN  products  could  explain  the  experimentally  observed  excitation  of  the 
HCN  bending  vibration  in  the  products  of  syw-triazine  i^otodissociated  at  193  nm  (Goates,  Chu,  and 
Flynn  1984).  The  lack  of  excitation  in  the  corresponding  C-H  and  C-N  bonds  could  be  explained  by  the 
small  changes  in  these  geometric  parameters  along  the  reaction  path  fiom  the  transition  state  to  the  HCN 
products. 

Since  such  concerted  triple  association/decomposition  reactions  are  considered  to  be  uncommon,  we 
attempted  an  analysis  that  could  show  why  these  reactions  occur,  using  projections  of  the  vibrational 
eigenvectors  of  syw-triazine  and  the  (HCN)3  cluster  onto  the  eigenvectors  at  various  points  along  the 
reaction  path  that  are  associated  with  the  direction  of  reaction  coordinate  for  the  concerted  reactions.  For 
points  all  along  the  reaction  path,  including  the  transition  state,  three  vibrational  modes  of  syw-triazine 
and  two  vibrational  modes  of  the  (HCN)3  cluster  project  strongly  onto  the  reaction  path.  Two  of  the 
modes  that  project  strcmgly  in  both  syw-triazine  and  the  (HO*l)3  cluster  are  symmetric  ring-breathing 
motions,  while  the  third  in  syw-triazine  is  a  C-H  rock.  Two  other  vibrations  of  the  (HCN)3  cluster  and 
a  vibration  in  syw-triazine  also  couple  to  portions  of  the  reaction  path,  though  not  as  strongly  as  the 
aforementioned  breathing  modes.  These  projections  indicate  that  certain  vibrational  modes  of  syw-triazine 
and  the  (HCN)3  cluster  are  coupled  to  the  reaction  coordinate,  providing  mechanisms  through  which 
reaction  energy  of  the  two  species  can  transfer  into  the  reaction  path,  resulting  in  the  ccmcetted 
association/decomposition  reactions.  The  existence  of  the  cyclic  (HCN)3  cluster  in  the  reaction  coordinate 
reduces  the  gross  entropic  effects  that  would  hinder  the  concerted  triple  association  reaction,  putting  the 
system  in  a  {uvferential  arrangement  for  concerted  association  if  reaction  energy  is  iqrpropriately 
distributed  among  this  cluster. 
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1  UNIVERSITY  OF  TEXAS 

DEPT  OF  CHEMISTRY 
ATTN  W  GARDINER 
AUSTIN  TX  78712 

1  VA  POLYTECH  INST  AND  STATE  UNIV 

ATTN  J  A  SCHETZ 
BLACKSBURG  VA  24061 
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NO.  OF 
COPIES 

1 

2 

1 

1 

1 

1 

1 

1 

1 

2 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


APPLIED  COMBUSTION  TECHNOLOGY  INC 
ATTN  A  M  VARNEY 
P  O  BOX  607885 
ORLANDO  FL  32860 

APPLIED  MECHANICS  REVIEWS 
ASME 

ATTN  R  E  WHITE  &  A  B  WENZEL 
345  E  47TH  STREET 
NEW  YORK  NY  10017 

TEXTRON  DEFENSE  SYSTEMS 
ATTN  A  PATRICK 
2385  REVERE  BEACH  PARKWAY 
EVERETT  MA  02149-5900 

BATTELLE 

TWSTIAC 

505  KING  AVENUE 

COLUMBUS  OH  43201-2693 

COHEN  PROFESSIONAL  SERVICES 
ATTN  N  S  COHEN 
141  CHANNING  STREET 
REDLANDS  CA  92373 

EXXON  RESEARCH  &  ENG  CO 
ATTN  A  DEAN 
ROUTE  22E 
ANNANDALE  NJ  08801 

GENERAL  APPLIED  SCIENCE  LABS  INC 
77  RAYNOR  AVENUE 
RONKONKAMA  NY  11779-6649 

GENERAL  ELECTRIC  ORDNANCE  SYSTEMS 
ATTN  JMANDZY 
100  PLASTICS  AVENUE 
PITTSFIELD  MA  01203 

GENERAL  MOTORS  RSCH  LABS 
PHYSICAL  CHEMISTRY  DEPARTMENT 
ATTN  T  SLOANE 
WARREN  MI  48090-9055 

HERCULES  INC 
ATTN  W  B  WALKUP 
E A  YOUNT 
P  O  BOX  210 

ROCKET  CENTER  WV  26726 


1  HERCULES  INC 

ATTN  R  V  CARTWRIGHT 
100  HOWARD  BLVD 
KENVIL  NJ  07847 

1  ALIiANT  TECHSYSTEMS  INC 

MARINE  SYSTEMS  GROUP 
ATTN  D  E  BRODEN  MS  MN50  2000 
600  2ND  STREET  NE 
HOPKINS  MN  55343 

1  ALLIANT  TECHSYSTEMS  INC 

ATTN  R  E  TOMPKINS 
MN  11  2720 

600  SECOND  ST  NORTH 
HOPKINS  MN  55343 

1  IBM  CORPORATION 

RESEARCH  DIVISION 
ATTN  A  C  TAM 
5600  COTTLE  ROAD 
SAN  JOSE  CA  95193 

1  nr  RESEARCH  INSTITUTE 

ATTN  R  F  REMALY 
10  WEST  35TH  STREET 
CHICAGO  IL  60616 

1  LOCKHEED  MISSILES  &  SPACE  CO 

ATTN  GEORGE  LO 
3251  HANOVER  STREET 
DEPT  52  35  B204  2 
PALO  ALTO  CA  94304 

1  OUN  ORDNANCE 

ATTN  V  MCDONALD  LIBRARY 
P  O  BOX  222 
ST  MARKS  FL  32355-0222 

1  PAUL  GOUGH  ASSOCIATES  INC 

ATTN  P  S  GOUGH 
1048  SOUTH  STREET 
PORTSMOUTH  NH  03801-5423 

1  HUGHES  AIRCRAFT  COMPANY 

ATTN  T  E  WARD 
PO  BOX  11337 
TUCSON  AZ  85734-1337 
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3 

1 

1 

1 

3 


3 

1 
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ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


SCIENCE  APPLICATIONS  INC 
ATTN  R  B  EDELMAN 
23146  CUMORAH  CREST 
WOODLAND  fflLLS  CA  91364 

SRI  INTERNATIONAL 
ATTN  G  SMITH 
D  CROSLEY 
DCiOLDEN 

333  RAVENSWOOD  AVENUE 
MENLO  PARK  CA  94025 

STEVENS  INSTITUTE  OF  TECH 
DAVIDSON  LABORATORY 
ATTN  R  MCALEVY  m 
HOBOKEN  NJ  07030 

SVERDRUP  TECHNOLOGY  INC 
LERC  GROUP 

ATTN  R  J  LOCKE  MS  SVR  2 
2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 

SVERDRUP  TECHNOLOGY  INC 
ATTN  J  DEUR 

2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 

THIOKOL  CORPORATION 
ELKTON  DIVISION 
ATTN  R  BIDDLE 
R  WILLER 
TECH  LIB 
P  O  BOX  241 
ELKTON  MD  21921 

THIOKOL  CORPORATION 
WASATCH  DIVISION 
ATTN  S  J  BENNETT 
P  O  BOX  524 

BRIGHAM  CITY  UT  84302 

UNITED  TECHNOLOGIES  RSCH  CENTER 
ATTN  A  C  ECKBRETH 
EAST  HARTFORD  CT  06108 

UNITED  TECHNOLOGIES  CORP 
CHEMICAL  SYSTEMS  DIVISION 
ATTN  R  R  MILLER 
P  O  BOX  49028 
SAN  JOSE  CA  95161-9028 


1  UNIVERSAL  PROPULSION  COMPANY 

ATTN  H  J  MCSPADDEN 
25401  NORTH  CENTRAL  AVENUE 
PHOENK  AZ  85027-7837 

1  VERTTAY  TECHNOLOGY  INC 

ATTN  E  B  FISHER 
4845  MILLERSPORT  HIGHWAY 
EAST  AMHERST  NY  14051-0305 

1  FREEDMAN  ASSOCIATES 

ATTN  E  FREEDMAN 
2411  DIANA  ROAD 
BALTIMORE  MD  21209-1525 

3  ALLIANT  TECHSYSTEMS 

ATTN  C  CANDLAND 
L  OSGOOD 
R  BECKER 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

1  DIRECTOR 

US  ARMY  BENET  LABS 
ATTN  AMSTA  AR  CCB  T 
SAM  SOPOK 
WATERVLIET  NY  12189 
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36 


ORGANIZATION 


ABERDEEN  PROVING  GROUND 


DIR  USARL 

ATTN:  AMSRL-WT-P,  A  HORST 

AMSRL-WT-PC, 

R  AFIFER 
GF  ADAMS 
W  R  ANDERSON 
R  A  BEYER 
S  WBUNTE 
C  F  CHABALOWSKI 
K  P  MCNEILL-BOONSTOPPEL 
A  COHEN 
RCUMPTON 
R  DANIEL 
D  DEVYNCK 
NFFELL 
BE  PORCH 
JMHEIMERL 
A  J  KOTLAR 
MRMANAA 
W  F  MCBRATNEY 
K  L  MCNESBY 
S  VMEDLIN 
MS  MILLER 
A  W  MIZIOLEk 
S  H  MODIANO 
J  B  MORRIS 
J  E  NEWBERRY 
S  A  NEWTON 
R  A  PESCE-RODRIGUEZ 
BMRICE 
R  C  SAUSA 
M  A  SCHROEDER 
J  A  VANDERHOFF 
MWENSING 
A  WHREN 
J  M  WIDDER 
C  WILLIAMSON 
AMSRL-CI-CA,  R  PATEL 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number/ Author  ARL-TR-1096  (Pail _ Date  of  Report  June  1996 _ 

2.  Date  Report  Received  _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.)  _ _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


DEPARTMENT  OF  THE  ARMY 


OFFICIAL  BUSINESS 

BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  0(I01,APQ,MO 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


DIRECTOR 

U.S.  ARMY  RESEARCH  LABORATORY 
ATTN:  AMSRL-WT-PC 

ABERDEEN  PROVING  GROUND,  MD  21005-5066 


